function a = test()
% TEST - 
%  
%    A=gausshermi(@cos,1,2,10)
%    B=sqrt(pi)*exp(-1/4)
%    a=A-B
%
    a = adp(@g,0,1)
    a-pi
end

function f = fun(x)
% FUN - 
%   

    f=cos(x);
end


function f = g(x)
% FUN - 
%   
    sx = size(x);
    f = ones(sx);
    
    for i = 1:sx(1)
        for j = 1:sx(2)
            f(i,j)=4/(1+x(i,j)^2);
        end
    end
    f=f;
end
